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ABSTRACT 

We present a stndy of the high rotational bands (J > 2) of H 2 toward 4 early 
type galactic stars: HD 73882, HD 192639, HD 206267, and HD 207538. In each 
case, the velocity dispersion - characterized by the spectrnm htting parameter b 
- increases with the level of excitation, a phenomenon that has previously been 
detected by the Copernicus and IMAPS observatories. In particular, we show 
with 4(7 conhdence that for HD 192639 it is not possible to £t all J levels with 
a single b value, and that higher b values are needed for the higher levels. The 
amplitude of the line broadening, which can be as high as 10 kms“^, makes ex¬ 
planations such as inhomogeneous spatial distribution unlikely. We investigate 
a mechanism in which the broadening is due to the molecules that are rotation- 
ally excited through the excess energy acquired after their formation on a grain 
(H 2 -formation pumping). We show that different dispersions would be a natu¬ 
ral consequence of this mechanism. We note however that such process would 
require a formation rate 10 times higher then what was inferred from other obser¬ 
vations. In view of the difficulty to account for the velocity dispersion as thermal 
broadening (T would be around 10 000 K), we conclude then that we are most 
certainly observing some highly turbulent warm layer associated with the cold 
diffuse cloud. Embedded in a magnetic held, it could be responsible for the high 
quantities of CH+ measured in the cold neutral medium. 
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1. Introduction 

The population distribution of the different rotational levels (J) of H 2 provides detailed 
information about diffuse and translucent clouds. The kinetic temperature, derived from the 
column density distribution of H 2 in the J = 0, 1 and 2 levels, is usually about 80 K. On 
the other hand, the excitation temperature obtained from the J = 3 - 6 levels is typically 
of several hundred K. The measurement of these column densities is an important probe 
of the physical conditions of the interstellar gas. Temperature, density, UV radiation held 
and other parameters can be inferred from such information, but a reliable model for the 
distribution of populations is required. 

There are three distinct mechanisms that likely determine the population of the excited 
levels: ultraviolet pumping, H 2 formation on grains, and high temperature collisional pro¬ 
cesses. Of these, the UV photoexcitation process (as described in Black & Dalgarno 1976) 
has been generally considered to be dominant. In the case of ( Oph (Black & Dalgarno 
1977), for instance, good matches were obtained for both the H 2 population distribution and 
the abundances of all known chemical species, except CH^. 

Additional information, explored in this paper, can be obtained from the presence of a 
measurable velocity dispersion which is an increasing function of the rotational energy level 
(J). This effect was first seen in several Copernicus observations, and was seen in both the 
curve of growth 6-value and in the line widths (Spitzer & Cochran 1973; Spitzer et ah 1974). 
More recently, Jenkins & Peimbert (1997) using IMAPS data (with a resolution power of 
120,000), showed, for the two main resolved components in the line of sight towards C, Ori A, 
a clear broadening of the H 2 lines increasing with rotational level. Along that line of sight, 
at least in the component showing the highest broadening, there is no doubt that the effect 
is linked to the excitation process of the H 2 molecule. 

It is this interplay between rotational excitation and velocity dispersion that we explore 
in this paper. Data from the FUSE FUV satellite (Moos et al. 2000) offers new insights 
on this issue. The wavelength range covered by FUSE (905 to 1187 A) contains more 
than twenty Lyman vibrational transitions, as well as six Werner bands, providing a large 
number of transitions for each rotational ground state and allowing measurements over a wide 
range of oscillator strengths. In addition, the high sensitivity of the instrument puts within 
observational reach some interesting lines of sight with high extinction, such as HD 73882, 
with E{B — V) = 0.73 (Snow et al. 2000). Although the spectral resolution does not allow 
us to resolve the different absorption components, we will show that saturation effects can 
also reveal the velocity dispersion of the lines. 

This paper is organized as follows. The observations and data reduction are described 
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in Section 2. Section 3 presents the H 2 analysis along the four sightlines listed in Table 1. 
Calculations were done using the curve of growth method (hereafter COG), and prohle htting 
(hereafter PF). Both methods are explained and discussed. In Section 4, we compare the 
observations with the theoretical explanation of the excitation, and discuss the consequences 
on the chemistry of the cloud. In the Appendix we describe a model of H 2 formation and 
subsequent cooling, which could explain the observed velocity dispersion. 


2. Observations &; Data Reduction 

The FUSE mission, its planning, and its in-orbit performance are discussed by Moos 
et al. (2000) and Sahnow et ah (2000). Briefly, the FUSE observatory consists of four 
coaligned prime-focus telescopes (two SiC and two LiF) and Rowland-circle spectrographs. 
The SiC gratings provide reflectivity over the range 905-1105 A, while the LiF have sen¬ 
sitivity in the 990-1187 A range. Each detector is composed of two micro-channel plates, 
therefore a gap of ~ 5 A divides each of our spectra into two pieces. The list of the four 
targets studied in this work and the log of the observations are presented in Tables 1 and 2 
respectively. All data were obtained with the source centered on the 30” x 30” (LWRS) 
aperture with total exposure times ranging from 4.8 ks (HD 192639) to 25.5 ks (HD 73882). 
All our datasets have a S/N ratio per pixel around 10. The data were processed with version 
2.0.4 of the CalFUSE pipeline. Corrections for detector background, Doppler shift, geomet¬ 
rical distortion, astigmatism, dead pixels, and walk^ were applied, but no correction was 
made for the fixed-pattern noise. The 1-D spectrum was extracted from the 2-D spectrum 
using optimal extraction^ (Horne 1986; Robertson 1986). Instead of co-adding the different 
segments of the spectrum, we used only the segments that appear to have the best correc¬ 
tion of the distortion effects. Those are SiC2A (930-990 A), LiFlA (990-1080 A), SiClA 
(1080-1088 A), and LiF2A (1090-1187 A). Below 930 A, the high reddening of our targets, 
removing most of the flux, does not allow for reliable measurements. 

After binning the data by 4 pixels (~ 7 kms“^), the processed data have a S/N ratio 
of nearly 20 per bin, and a nominal spectral resolution of ~ 20 kms“^ (FWHM). 


^The FUSE detector electronics happens to miscalculate the X location of photon events with low pulse 
heights. This effect is called “walk”. 


^http: / / fuse.pha.jhu.edu/analysis/lacour/ 
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3. H 2 Measurements 

The FUSE wavelength range allows us to access a large number of H 2 absorption lines, 
corresponding to a wide range of rotational excitations. For each of the levels that we focus 
on (J = 2 — 7), we have measured, when available, column densities and b values. To ensure 
consistency of the measurements, we used two different methods to determine N and b, 
described below. 


3.1. Curve of Growth Method (COG) 

The measured equivalent widths (EqW) of each line studied in this work are summarized 
in Table 3. The stellar continuum in the vicinity of each line was estimated using a low-order 
Legendre polynomial ht to the data. The la error bars were computed taking into account 
four types of errors, added in quadrature: 

• The statistical errors, supposed to be a white Poissonian noise. These errors (roughly 
the square root of the count rate) are computed by the pipeline for each pixel. The 
total statistical error over each line is therefore the quadratic sum of the error of each 
integrated pixel (more information can be found in Appendix A of Sembach & Savage 
1992). 

• The background uncertainties, proportional to the exposure time. They have been 
estimated by the FUSE science data processing team to be at the level of 10 % of 
the computed background^. This error is calculated by the pipeline, and added to the 
statistical errors. 

• The continuum placement error, which depends mainly on the S/N ratio in the vicinity 
of the line. To estimate this error, we shift the continuum by 1 to 3 % (depending 
on the S/N ratio), calculating a lower and upper value for the EqW. The difference is 
taken as the la error. We note here that because of to the stellar type of the targets 
(see Sp.T. in Table 1), very few stellar lines are present, and are easily distinguishable 
with the interstellar line due to their thermal broadening. 

• The systematic uncertainties, which are the most difficult errors to quantify. They may 
come from geometrical distortions, walk, dead pixels, point spread function (PSF), hxed 
pattern noise, etc. Most of these distortions are corrected by the pipeline, but these 


^ htt p: / / fuse. pha. j hu. edu / analysis / calfuse_wp3. html 
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effects may nevertheless have a non-negligible influence on our measurements. More¬ 
over, there is no way to estimate the effect over a single absorption line. Assuming that 
systematic errors are homogeneous over our measurements, we adjusted the systematic 
errors to be proportional to the EqW. The factor of proportionality is set so that the 
total of the COG £t is equal to the number of degrees of freedom. To avoid any 
bias, a proportionality factor was obtained independently for each sight line and each 
species (i.e., for each J level), which is possible because the number of spectral lines 
being measured is statistically significant. The resulting factors are in the range of 1 
to 8 %. 

To determine column densities and b values, we fit a set of 300 single-Gaussian curves 
of growth to our measured EqWs. They were obtained by integrating a Voigt profile over a 
large number of b values and damping factors. For each species, was calculated for each 
b value and column density {N). 

The best COG hts are shown in Figure 1. The upper plots show the resulting COG 
using a single b value for all the rotational levels. The value of compared to the number of 
degrees of freedom is the best mathematical tool to evaluate the goodness of a £t. For these 
hts, the values are 52/35 for HD 73882, 110/58 for HD 192639, 54/44 for HD206267, and 
37/36 for HD 207538. The probabilities of having xt equal to or larger than those values are 
3% (> 2a), 0.0045% (> 4(j), 14% (> la) and 42% respectively. Therefore, the result is not 
signihcant for HD 207538 but HD 192639, and to a lesser extent, HD 73882 and HD 206267, 
clearly have an inconsistency in the £t. The middle plots display our best hts using diherent 
b values for each rotational level and the bottom plots show the Ay^ as a function of b, 
again, for each rotational level. All errors listed here are at the 2a level, corresponding to a 
Ay^ of 4. 

To check the possibility that some systematic problem in one (or a few) lines could 
induce a “false broadening ehect” on the HD 192639 dataset, we randomly removed half of 
our measurements, then the other half. The results, summarized in Table 4, conhrm the 
presence of the broadening. 


3.2. Profile Fitting Method 

The spectral resolution (si 20 kms“^, varying over several kms“^, depending on the 
detector used), is insufficient to directly determine variations in the broadening of lines, 
which are typically of the order of a few kms“^. Instead, it is the relative shape of the lines, 
combined with knowledge of the oscillator strengths, that provides meaningful information. 
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We used a profile fitting routine, Owens (Lemoine et al. 2002; Hebrard et al. 2002), developed 
by M. Lemoine at the Institut d’Astrophysique de Paris, which allows htting all the lines 
simultaneously while considering different b values for each rotational level. To minimize 
systematic errors induced by the PSF, we allowed the width of the PSF to vary by ~ 25 % 
about the nominal value (15-25 kms“^). An important advantage of htting multiple species 
is the added ability to work with semi-blended lines. The lines that were included with this 
method are listed in Table 3. This was particularly helpful in constraining the J = 2 b value, 
because of the weak oscillator strength of the 1112.5 A blended line. Figure 2 shows some 
sample best hts for the 1017.8 A (J = 5), 1116.0 A (J = 4) and 1115.91 A (J = 3) lines. 
In Figure 3 we plot the for the b values used in the hts. To account for systematic 
uncertainties, we scaled the errors (by a factor of 1 to 2) so that the minimum is equal 
to the number of degrees of freedom (for more on the y^ technique while using Owens, see 
Hebrard et al. 2002). 

To convince ourselves that prohle htting with a single b value for all levels does not 
accurately describe our data, we did the htting for HD 73882 and HD 192639 and compared 
the htting on one particular line. The 2 panels in the right column of Figure 4 show the best 
ht for the 1017.8 A line (segment LiFlA) when a single b value was used for all rotational 
levels, while the left panels show the best hts using individual b values. Without taking 
systematic uncertainties into account (i.e., without scaling the errors), the y^ for the left 
panels are close to the number of degrees of freedom (28/27 for HD 73882 and 24/21 for 
HD 192639), indicating that the hts are reliable. On the other hand the y^ for the right 
panels (62/27 and 44/21 respectively) indicate an inconsistency. Note that the PSF was a 
free parameter and had a value of 2.96 binned pixels (23.4 kms“^) for HD73882 and 2.98 
binned pixels (23.2 kms“^) for HD 192639 (left column). For the right column, the obtained 
PSFs are 2.67 pixels (21.1 kms“^) for HD 73882 and 2.42 pixels (19.1 kms“^) for HD 192639. 
The diherences in the PSFs are within the FUSE resolution uncertainties (see Moos et al. 
2000 ). 


3.3. Results 

Tables 5 list the column densities and the b values, respectively, derived using the 
methods described above (2(T uncertainties). Prohle htting allows us to quickly determine 
upper limits on column densities for J levels which only have transitions that are too weak 
to be used with the COG method. Column densities for the J = 0 and 1 levels, for the 
four stars, are from Rachford et al. (2002). The b values are consistent between PF and 
COG. In each case, but with diherent reliability levels, the velocity dispersion shows the 
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same increasing trend with the H 2 excitation levels. Since b values and column densities are 
interdependent, it is an important result for observers investigating saturated H 2 lines. If 
the broadening of the higher J levels is assumed to be equal to those of the lower J levels, 
then there is the possibility of considerably overestimating the column densities. 


4. On the source of excitation and broadening 

In light of the results presented above a question must then be asked: What causes the 
increase in the broadening of the H 2 lines with increasing J level? An obvious explanation 
would be that we may have a broad component revealed at high rotational level by an high 
excitation temperature. To test this possibility, we plot in Figure 5 the curves of growth for 
a sightline having two components with different velocity dispersions. Because the relative 
strength differs from one rotational level to the other, each level correspond to a curve of 
growth with a different shape. As an example, we used the excitation diagram of target 
HD 193639 to £t two components which are associated with the J = 2 level for one, and the 
J = 5 level for the other. We then fitted the EqWs on the COGs of their corresponding 
excitation level. The fit explain why an effect due to the variation in the ratio between two 
components, one broad and the other narrower, can be seen as a variation in the broadening. 
We note nevertheless that this explanation is incompatible with Component 1 observed with 
IMAPS towards C, Ori A (Jenkins & Peimbert 1997). In this case, the third rotational level 
column density is so low, that it cannot belong to a different component. 

Finding an explanation for the presence of a broad component is another challenge. It 
may be the key behind the source of both the excitation of H 2 and the presence of large 
amounts of CH'*". Specifically, the fast ion-molecule reaction, CH’*' J- H 2 — CH^ J- H, in 
cold gas, predicts CH'*' column densities far below the observed levels (Watson 1974). This 
is also the case toward our sightlines (A^(CH''') > 10^^ cm“^). The solution might be in 
warm interstellar gas (T > 10^ K) in which the endothermic reaction C’*' J- H 2 —> 

H - 0.4eV can provide an equilibrium density close to the observed levels. The presence 
of a warm component received strong support by the observation of a correlation between 
CH'*' and the rotationally excited H 2 (Lambert & Banks 1986). However, since until now 
no direct observations of this warm component have been obtained, parameters such as its 
density and temperature are unknown. The increase of b with increasing J level seems to 
be direct evidence of a warm component. Several excitation mechanisms, discussed below, 
could be responsible for this effect. 



4.1. UV-pumping 


The ro-vibrational cascading releases its energy through infrared photons (Black & 
Dalgarno 1976). Like photoexcitation, such energy loss does not change the kinetic energy 
of the molecules, and therefore it does not affect the velocity dispersion. Heating of the gas 
can nevertheless occur, through photodissociation of H 2 and photoelectron emission from 
dust grains. However, such process would require a high UV held (as the one towards the 
Pleiades cluster, e.g. White 1984), and is an unlikely explanation for a broadening of up to 
10 kms“^. 


4.2. H 2 -formation pumping 

When molecules are created on the surface of grains, they carry away most of the initial 
energy (?« 4.5 eV) which provides the kinetic, rotational and vibrational excitation of H 2 . 
Support for this mechanism was obtained by Wagenblast (1992) who calculated the column 
density ratio between the J = 4 to 7 levels in good agreement with observations. However, 
absolute column density calculations do not exist. To address this, we constructed a time- 
dependent model in which we followed the stochastic evolution of the molecules after their 
grain formation. Details of the model are given in the Appendix. According to this model, 
dispersions appear as a natural consequence of the equilibrium among the various excitation 
and decay processes. Figure 6 shows the expected broadening, as a function of the density 
and the rotational level. Using the expected broadening in conjunction with Equation A3 we 
are able to calculate the densities and formation rates needed to explain both the velocity 
dispersions and the column densities of the broader rotational levels (J > 4). The densities 
are in agreement with previous analysis of the C I fine-structure excitation for HD 192639 
and HD 206267, which led to to estimated densities of 16 cm“^ (Sonnentrucker et ah 2002) 
and 30 cm“^ (Jenkins & Tripp 2001), respectively. However, the formation rates implied by 
our models (i? in Table 6), do not agree with previously calculated values {R ^ 3 x 10“^^ 
cm^s“^ in Jura 1975; Gry et ah 2002). There is a factor of approximately ten between the 
rate needed to explain the column densities and the rates mentioned above. 

A second argument against H 2 -formation pumping as responsible for the broadening 
of the excited states comes from the difficulty to account for the CH^ column density. We 
report on the upper right panel of Figure 7 the threshold energy for the C^-|-H 2 reaction, and 
found that the time during which the molecule is kinetically warm is fwarm = 4 x 10®/n s (the 
cooling time is roughly inversely proportional to the density n, in cm“^). Assuming that the 
density is spatially uniform along the sightline, we can obtain the column density of warm H 2 
as a function of the atomic hydrogen and the formation rate : A^(H 2 )warm = A^(H) ni? fwarm- 
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Considering the formation rates from Table 6, we obtain a ratio iV(H 2 )warm/-^(H) ~ 1.2 x 
10“^, a value several orders of magnitude below what is needed to explain the column density 
of CH+ (iV(H 2 )warm/^(H 2 )cooi ~ 10“^ in Lambert & Banks 1986). 

4.3. Collisional excitation in a warm environment 

Warm low density interstellar gas surely is present along the lines of sight. Field et al. 
(1969) were the hrst to show that warm gas could be thermally stable at low densities. Such 
gas (uh ~ 0.1 cm“^; T ^ 7000 K) appears to be the major constituent - in mass - of 
the local interstellar cloud (Frisch & York 1983; Linsky 1996; Lallement 1998), but a recent 
survey of the LISM, performed by Lehner et al. (2003), showed that in this medium, the H 2 
molecular fraction is low, close to 10“^. The problem appears to be that at such densities 
and temperatures, H 2 formation on grains becomes negligible. Other routes exist, such as 
formation through hydrogen ions (Black et al. 1981), but the rates for the process are low, 
as is the observed amount of such ions (Andre et al. 2002). The same arguments can be 
used against the presence of H 2 in the warm postshock gas of dissociative shocks. However, 
Jenkins & Peimbert (1997) concluded anyway that in the line of sight towards ( Ori A, the 
dispersion of the J = 5 line could be the trace of an ongoing J-type shock after which the 
gas recombines and cools. It is true that this theory was corroborated by a shift of the 
line centers - the higher J levels being shifted towards lower velocities - which is difficult to 
explain otherwise. 

But shifts could also be the trace of slower shocks, C-type shocks. Even though the 
conditions in the postshock gas are not favorable in terms of H 2 formation, the fact that 
it is non-dissociative makes it possible to contain a large amount of heated H 2 . Hence, 
Elitzur & Watson (1978) showed that such shocks could heat a large portion of the cloud 
with temperatures of several thousands Kelvin. This mechanism would generate both the 
broadening and the CH+ column densities. However, there are several arguments against 
the presence of a single important shock front: predictions of Y(OH) produced via an 
endothermic reaction would be signihcantly higher than what is observed, and signihcant 
velocity shifts required for this type of shocks are usually undetected. We can however derive 
the environment parameters in the hypothesis of a thermal broadening of the postshock gas. 
We listed in Table 7 the kinetic temperature (under the assumption that the linewidths 
are mostly due to thermal broadening), and the excitation temperature (derived from the 
J = 5 and 7 levels) for the four lines of sight. The difference between the two temperatures 
constrain the ratio between collisional excitation and radiative de-excitation rates. Hence, 
we used the rates from Le Bourlot et al. (1999) and Wolniewicz et al. (1998) to infer the 
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densities, and therefore the pressure log{P/k) = log{nT). We note that we are not able to 
derive lower limits of the J = 7 rotational column density. However, we used the best htting 
values to infer the pressure of the neutral gas, and obtained values signihcantly higher than 
those derived from the C I hne structure; log(P/k) ^ 3.1 for HD 192639 (Sonnentrucker et al. 
2002) and 3.5 for HD206267 (Jenkins & Tripp 2001). 

From the pressure calculations, we argue that we are more likely in the presence of 
material cooler than what can be derived from a thermal velocity dispersion. The molecule 
must still be hot enough to explain the excitation temperature (approximatively a thousand 
Kelvin), but not up to the temperature required for a thermal broadening (T^ in Table 7). 
Hence, the only remaining explanation of the velocity dispersion is that we are looking at a 
warm and turbulent layer, most probably intimately associated with the cold medium. To 
explain the presence of such layer, some invoke supra-thermal velocities of ions relative to 
the neutrals driven by multiple magneto-hydrodynamic (MHD) criss-crossing shocks (Gredel 
et al. 2002). Others invoke the intermittency of turbulence and the existence of localized tiny 
warm regions, transiently heated by bursts of ion-neutral friction and viscous dissipation in 
coherent and intense small small vortices (Joulain et al. 1998). The main physical difference 
between each phenomenon is the thickness and the crossing time. While a warm MHD 
postshocks layer can have a thickness of ~ 0.1 pc, coherent vortices threaded by magnetic 
helds may have radii as small as ~ 20 AU. Both would achieve peak temperature around 
lOOOK, but with differential velocities around 10 kms“^ for the MHD shocks and around 4 
kms“^ in a coherent vortice. 

Towards our targets, the situation would be perfectly describe by both explanations. 
The many small-scale shocks or vortices would create the observed velocity dispersion, while 
embedded magnetic helds would generate differential velocities between ion and neutral 
species, reacting into CH"*" through the C’*' J- H 2 —> CH+-I- H endothermic reaction. It 
would also give an answer to multiple quests for a warm layer (e.g. CO, HCO-I-, H 2 O in 
Pety & Falgarone 2000; Liszt & Lucas 2000; Neufeld et al. 2002). However, other heating 
process may be possible, and are difficult to rule out due to our low spectral resolution. 
The probability distribution functions of these high rotational levels — either in FUV or IR 
(Verstraete et al. 1999, Falgarone et al. in press) — would eventually help us to distinguish 
between one process and the other. 


5. Summary 

We observed four highly reddened {E(^b-v) ~ 0.6) lines of sight in which the saturation 
of the higher rotational levels of H 2 allows us to infer their velocity dispersion. We measured 
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broadenings up to 10 kms“^, increasing with the energy of the rotational level. Considering 
the fact that it was already observed towards several other sightlines (Spitzer & Cochran 
1973; Spitzer et ah 1974; Jenkins & Peimbert 1996), we suggest that this phenomenon is a 
fairly common one. As a hrst result, we suggest caution in investigating saturated H 2 lines 
since assuming an identical linewidth for all the rotational levels could lead to signihcant 
systematic errors on the column densities. 

We looked at the possible sources of rotational excitation to see which could induce 
such broadening of the absorption lines. We ruled out UV-pumping which can not create 
a velocity dispersion by itself. We constructed a time dependent model of the state of the 
molecule following formation on the grain. From it, we deduced that such mechanism may 
be responsible for the broadening, but would however need a formation rate ten times the 
one derived from previous studies (Jura 1975). We also note that the amount of warm H 2 
created would not be enough to account for the observed column densities of CH'*'. 

We conclude that the more likely explanation is the presence of a turbulent, warm layer 
in the molecular cloud. The temperature would need to be over 600 K to account for the 
rotationally excited H 2 , and the velocity dispersion of the gas patches around 8 kms“^. 
Small criss-crossing shocks or vortices could be the phenomenon behind this turbulent layer. 
Magneto-hydrodynamic waves created inside them would convert the kinetic energy to create 
the observed amount of CH'^ through its endothermic reaction with H 2 . 


The authors thank N. Balakrishnan for his work on the modeling of the cooling of hot 
H 2 . The authors also thank Ed Jenkins for invaluable discussions, thorough comments on 
an early draft and analysis of Cl I lines which do not hgure in this paper. Further thanks go 
to E. Falgarone, B. Racheford, P. Sonnentrucker and N. Lehner for helpful comments and 
discussions. SL would also like to personally thank the entire JHU FUSE Team for their 
hospitality, their knowledge and their kindness in sharing it. Part of this work has been done 
using the prohle htting procedure Owens.f developed by M. Lemoine. French participants 
are supported by CNES. This work is based on data obtained for the Guaranteed Time Team 
by the NASA-CNES-CSA FUSE mission operated by The Johns Hopkins University. 


A. Modeling H 2 formation 
A.l. Description 

In this section we describe our numerical model of H 2 excitation and de-excitation 
following H 2 formation on grains. We use it to show the existence of a possible correlation 
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between the excitation level of the molecule and the line broadening associated with that 
level. We also use it to calculate the amount of excited H 2 obtained from H 2 -formation 
pumping. Note that throughout most of this section we are dealing only with the gas that is 
excited by formation pumping. The effects described here are not affected by the presence 
of another gas component (at the same level) excited by a different mechanism. 

The formation of H 2 on grains has been discussed at length in the literature (e.g. Gould 
& Salpeter 1963; Knaap et ah 1966; Augason 1970; Hollenbach & Salpeter 1971; Lee 1972). 
Our model is based on the fact that when H 2 molecules are formed, they will leave the surface 
of the grain carrying with them excess energy, which will be distributed as surface, kinetic, 
rotational and vibrational energy. If we can determine the cross-sections for the various 
processes by which the molecules either lose or convert energy from one form into another, 
we can then estimate the average kinetic energy over the lifetimes of each species. Adding 
in the formation rates, we are able to predict the column density and velocity dispersion of 
each level. 


A. 1.1. Initial distribution 

Recent simulations of H 2 formation have been performed for graphite surfaces (Parneix 
& Brechignac 1998; Meijer et ah 2001) and icy interstellar dust (Takahashi et ah 1999). In 
all cases, the molecules leave the grain with an initial energy of 4.48 eV, but the distribution 
among binding, kinetic, and rotational energy differ considerably. We based our model on the 
quasi-classical computer simulation of the Parneix & Brechignac paper, due to their choice 
of computational method and physical model, which includes an entrance-channel barrier on 
the potential energy surface of the grain. The initial distribution over the vibrational levels 
corresponds to Fig. 11 of their paper, and averages around 1.1 eV, with a peak at n = 0. 
The distribution of rotational energy, with a mean value of 0.7 eV, corresponds to the Fig. 9 
from the same paper, and is the same for each vibrational level. 1.7 eV is carried as kinetic 
energy, and the remainder is used to desorb the molecule from the grain. The implications of 
these choices are discussed in the “Model results & considerations” Section. The ortho/para 
repartition (OPR) is an important factor in the column density calculations, but fortunately, 
it does not affect the results concerning the velocity dispersions. In view of the results of 
Persson & Jackson (1995), we decided not to include the statistical weight factor of 1/3. 
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A. 1.2. Kinetic cooling 

Starting with an initial kinetic energy, the H 2 molecule will cool down, while interacting 
with its environment. Because the formation rate is proportional to the density of atomic 
hydrogen, most of the formation of molecular hydrogen takes place in the photodissociation 
region, where hydrogen is mostly atomic. Moreover, simulations (Le Petit et al. 2002) show 
a sharp decrease of the molecular fraction as we go deeper in the cloud, implying a formation 
rate only marginal inside the cloud. It leaded us to approximate the formation environment 
as being purely atomic. Under this assumption, we needed just three parameters to compute 
the kinetic cooling: the H-H 2 cross-section (a), the mean energy loss per collision ( 7 ), and 
the density (n). While n is just a property of the medium, the two other parameters depend 
on the energy of the molecule. To obtain a, we extrapolated to higher energies the values in 
Fig. 2 of Clark & McCourt (1995). The results are plotted in the upper left panel of Figure 7. 
7 was generated using equations (12) and (25) in Kharchenko et al. (1998), while n is left as 
a free parameter, varying from 0.1 to 1000 cm“^. The right upper panel of Figure 7 displays 
the estimated kinetic energy as a function of time for a medium with density n = 10 cm“^. 
One can see that the translational energy is lost fairly slowly, allowing time for the molecule 
to undergo endothermic reactions. For example, if we consider the threshold energy for the 
-|- H 2 —CH’''-|- H - 0.4eV reaction (dashed line), the molecule would have such or more 
energy during = 4 x 10® seconds (~ 10 ans). 


A. 1.3. Radiative cooling 

From the initial distributions, we generated a table of 20 rows (corresponding to the 
20 first rotational levels) and 7 columns (corresponding to the 7 first vibrational levels). In 
steps of 10® seconds, and using the radiative decay table from Wolniewicz et al. (1998), we 
computed the density of each excitation level from its formation until 3.5 x 10^^ seconds. 
The panel in the center of Figure 7 shows the densities (normalized to 1) of each vibrational 
level (independent of the rotational level). After a fairly short time (~ lO^s), all of the 
molecular hydrogen is in the ground vibrational state. Then, the model clearly shows the 
progression towards lower J levels as time increases (lower panels of Figure 7). The rising 
parts of the density curves are caused by the cascading down from higher J or n levels, while 
the steep drops occur at the radiative lifetime. Interestingly, these lifetimes (e.g., ~ 10^ 

s) are of the order of the kinetic cooling time (see previous paragraph). It follows that, 
depending on the density of the cloud, some rotational levels will be kinetically hot, while 
lower levels will not. 
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A.1.4. Inelastic collisional cooling & excitation 

In addition to the decay rates, collisional excitation and de-excitation may have an 
important effect, especially while the molecule is still highly energetic, or when the density 
is high. We used the rates from the web site http://ccp7.dur.ac.uk/cooling_by_h2/index.htnil 
(Forrey et ah 1997; Le Bourlot et ah 1999), in addition to the decay rates of the above section, 
to model the time dependent density distribution. The lower panels of Figure 7 show the 
resulting densities (normalized to 1) of the pure rotational levels. The specific conditions 
for this plot are a density of 10 particles per cm^, and an ambient temperature of 100 K. 
As discussed in the previous paragraph, it is apparent that the fraction of each excited state 
falls off rapidly at a time corresponding to the decay rate, and that inelastic collisions are 
not an important de-excitation mechanism. Numerical tests show that for densities below 
several hundred cm“^, inelastic collisions are negligible compared to radiative decay. 


A.2. Model results &; considerations 

Integrating the normalized column densities over time allows us to extract two pieces 
of information : 


1. The mean kinetic energy of each rotational levels, which can be directly linked to 
the velocity dispersion as long as we assume that the excitation is mainly due to its 
formation on the grain. The mean energy is obtained by: 

. _ Density(J,n,t) x E{t)dt 
Density {J,n,t)dt 



In practice the upper limit of the integration was set to 3.5 xlO^^ seconds, by which 
time the excited species have largely decayed. The velocity dispersion, i.e., the broad¬ 


ening, is therefore 



(A2) 


We plotted in Figure 6, the broadening of each rotational levels between 2 and 6, 
clearly showing that it is increasing with the rotational level. Again, this is only 
the case when the species are populated by H 2 -formation pumping. The broadening of 
higher rotational levels (J = 4 and 5) were htted on the curves, and estimated densities 
are reported for our four targets in Table 6. 


2 . Still assuming no other source of excited H 2 , we are also able to calculate the molecular 
hydrogen column densities of the excited levels. The densities depend on the formation 
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rate (ni?), and the atomic hydrogen density. Assuming that the density is spatially 
uniform along the line of sight, we can use the column density of atomic hydrogen 
A^(H), and therefore, have the following equation: 

r+oo 

A^(H 2 )(J, n) = A^(H) ni? / Density(J, n, (A3) 

Jo 

To remove the effect of the formation OPR on the column density, we considered the 
sum of the column densities of the J = 4 and 5 levels. The formation rates needed to 
create such column densities are reported in Table 6. 

The distribution of the formation-energy of H 2 has been calculated through many other 
theoretical models, yielding quite different results in all the different parameters. Two pa¬ 
rameters have a large influence on our results: the initial kinetic and rotational energies. 
On one hand, the kinetic energy will decide the amplitude of the velocity dispersion. For 
instance, whereas our model starts with a kinetic energy of 1.7 eV, Meijer et al. (2001) hnd 
an initial energy of 1.18 eV which would scale down all our broadening calculations by a 
factor of 1.2. More specihcally, from Equation A2, a broadening of 8 kms“^ can be explained 
if the kinetic initial energy is equal or above 0.95 eV. This is compatible with most but a 
few of the theoretical formation models. On the other hand, the initial rotational energy 
does constrain the calculated column density of the rotational levels. Because the radiative 
lifetimes of the higher levels are significantly smaller than that of the lower rotational lev¬ 
els, a first approximation can be that each level is populated by its initial population plus 
the initial population of the higher level. Therefore, any model yielding an average initial 
rotational energy above 0.4 eV would give similar column density results for J < 7. Except 
for Duley & Williams (1986), most of the theoretical models agree with a rotationally hot 
initial distribution, hence in accordance with our chosen model. 
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Fig. 1.— Curve of growth analysis for H 2 along the HD 73882, HD 192639, HD 206267 and 
HD 207538 lines of sight. The plots at the top assume a common 6-value for all J levels. In 
the middle plots, each J level is £t with a different 6-value. The plots in the bottom of the 
figure present as a function of 6 for the COGs displayed in the middle of the figure. In 
the plots, the different rotational levels are represented with dotted (J = 2), dashed 
(J = 3), dot-dashed (J = 4), and dot-dot-dot-dashed (J = 5) lines. Note that due to the 
non-linearity of the COG, the Ay^ sometimes shows two minima. In these cases, the second 
minimum is ruled out on the grounds that the resulting column densities are not physical. 
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Fig. 2.— Best profile fits to the J = 5 1017.8 A line (left) and J = 3 1115.9 A and J = 4 
1116.0 A lines (right), for the four sightlines. Each rotational level is £t with a different b 
value. The continuum is normalized to 1. 
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Fig. 3.— curves for the different b values used with profile fitting. The different 
rotational levels are represented with dotted (J = 2), dashed (J = 3), dot-dashed (J = 4), 
and dot-dot-dot-dashed (J = 5) lines. 
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HD 73882 



HD 192639 



HD 73882 



HD 192639 



Fig. 4.— Best profile fits of the J = 5 1017.8 A line for HD 73882 (top) and HD 192639 
(bottom) using a single b value for all J levels (right) and allowing different b values for each 
J level (left). The inadequacy of the single b value hts for all J levels is apparent both from 
close inspection of the hts and from the reduced values (see section 3.2). 



23 




Fig. 5.— Left panel : Excitation diagram of the H 2 along the HD 192639 line of sight. We 
htted the column densities with several components, with different excitation temperatures 
and velocity dispersions. The T^x = Tj= 2-4 component is supposed to have a broadening 
equal to the J = 2 species, and the T^x = T 7 = 5-7 component a broadening equal to the 
J = 5 species. Right panel: Two components curve of growth for J = 2 to 6 , with a relative 
strength in accordance with the two components htted on the higher levels of the excitation 
plot. A variation of the 6 -values appear due to the presence of the two components. 
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Fig. 6.— Line broadening (6) derived from the mean kinetic energy of each level following 
formation on a grain. Since the value of the kinetic energy is dominated by collisional 
cooling, the broadening is a strong function of the density (n). The thin vertical lines show 
the density towards HD 192639 (Sonnentrucker et ah 2002) and HD 206267 (Jenkins & Tripp 
2001 ). 
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Fig. 7.— Time dependent model of molecular hydrogen following formation on grains (see 
Appendix). Upper left panel: Effective H-H 2 cross-section as a function of the kinetic energy. 
Upper right panel : Kinetic energy of the H 2 molecule as a function of time. The dashed line 
corresponds to the energy threshold from which C’*' can react with H 2 to form CH’*'. Three 
lower panels : Distribution of H 2 molecules over the vibrational and the ortho and para 
rotational states, as a function of the molecule’s lifetime. Note that radiative vibrational 
transitions ensure decay towards the vibrational ground state on short time scales. The 
initial rise of the rotational density curves is due to the cascading of molecules down from 
the higher J and V levels, whereas the steep drops occur at the radiative lifetime of the 
level. Because collisional excitation and de-excitation are part of our simulation, all these 
plots depend on the density, here equal to 10 cm“^. 

















Table 1. 


Sightline and Stellar Properties 


Star 

1 (") 

bO 

V (mag) 

E{B - V) ^ 

Av^ 

Sp.T. 

HD 73882 

259.83 

+0.47 

7.27 

0.72 

2.28 

09HI 

HD 192639 

74.90 

+1.48 

7.11 

0.66 

1.87 

08e 

HD 206267 

98.98 

+3.71 

5.62 

0.52 

1.37 

06e 

HD 207538 

102.86 

+6.92 

7.30 

0.64 

1.43 

BOV 


’^Extinction parameters from the FUSE H 2 Survey (Rachford et ah 







Table 2. 


Log of FUSE Observations 


Star 

FUSE ID'’ 

Observation 

Date 

Number of 
Exposures 

Exposure Time 
(Ks) 

S/N ’’ 

HD 73882 

P1161301 

2000.01.24 

6 

11.9 

5.1 


P1161302 

2000.03.19 

8 

13.6 

4.6 

HD 192639 

P1162401 

2000.06.12 

2 

4.8 

8.1 

HD 206267 

P1162701 

2000.07.21 

3 

4.9 

10.2 

HD 207538 

P1162902 

1999.12.08 

4 

7.7 

6.2 


P1162903 

2000.07.21 

10 

11.2 

7.1 


Archival root name of target for FUSE PI team observations. 

’’Average per-pixel S/N for a 1 A region of the LIF la spectrum near 1070 A. 







Table 3. H 2 equivalent width measurements 


Species 

A (A) 

log (/A) 

HD 73882 WA(mA) 

HD 192639 WA(mA) 

HD 206267 WA(mA) 

HD 207538 WA(mA) 

H 2 j - 2 . 

941.606 

0.498 


111.4 ± 9.7 

88.3 ± 8.0 



957.660 

0.661 


148.1 ± 13.5 

115.1 ± 6.1 

147.6 ± 13.6 


975.351 

0.810 


151.1 ± 11.0 

124.9 ± 6.7 

144.7 ± 10.3 


1005.40 

0.998 


212.9 ± 11.5 

162.4 ± 7.7 

208.1 ± 10.8 


1016.47 

1.016 

323.6 ± 33.3 

245.6 ± 12.7 

183.2 ± 8.4 

229.1 ± 10.1 


1040.37 

1.030 

359.7 ± 31.9 

275.7 ± 13.5 

202.9 ± 8.8 

263.6 ± 13.5 


1053.29 

0.980 

283.0 ± 23.6 

273.5 ± 13.1 

193.7 ± 9.3 

270.1 ± 13.6 


1066.91 

0.881 

280.5 ± 23.9 

257.7 ± 11.6 

192.9 ± 9.2 

239.7 ± 10.4 


1081.27 

0.709 

256.5 ± 24.9 

200.1 ± 11.0 

162.8 ± 6.9 

202.6 ± 11.4 


1096.45 

0.417 

220.1 ± 27.4 

175.1 ± 8.5 

PF 

PF 


1112.51 

-0.111 

PF 

PF 

PF 

PF 

H 2 j - 3 . 

934.800 

0.820 


106.3 ± 9.3 




942.970 

0.729 


108.5 ± 14.2 

85.2 ± 8.5 



944.337 

0.519 


95.2 ± 13.9 

67.5 ± 9.4 



958.953 

0.930 


95.8 ± 8.6 

97.5 ± 6.5 

90.8 ± 7.7 


960.458 

0.674 


92.7 ± 8.5 

90.5 ± 6.8 

86.7 ± 5.9 


995.974 

1.218 


143.9 ± 9.6 

117.6 ± 8.3 

109.2 ± 6.6 


997.830 

0.942 


123.1 ± 11.9 

97.5 ± 9.4 

91.4 ± 8.9 


1006.42 

1.199 

PF 

PF 

PF 

PF 


1017.43 

1.270 

PF 

PF 

PF 

PF 


1019.51 

1.030 

195.0 ± 16.4 

116.0 ± 6.9 

113.1 ± 7.2 

101.8 ± 4.5 


1028.99 

1.250 


144.9 ± 10.4 


128.2 ± 5.3 


1041.16 

1.216 

204.0 ± 12.5 

159.3 ± 9.5 

146.3 ± 7.7 

122.9 ± 6.6 


1043.51 

1.052 

182.1 ± 11.1 

139.5 ± 8.1 

131.1 ± 6.6 

117.6 ± 4.8 


1053.98 

1.150 

216.0 ± 13.4 

159.5 ± 10.6 

137.6 ± 7.7 

131.7 ± 6.7 


1056.48 

1.006 

192.8 ± 12.0 

125.6 ± 6.7 

122.6 ± 7.5 

121.6 ± 5.5 


1067.48 

1.028 

176.0 ± 15.7 

129.9 ± 7.6 

130.9 ± 7.4 

119.0 ± 5.3 


1070.15 

0.909 

196.0 ± 12.5 

134.6 ± 8.1 

129.0 ± 6.9 

110.0 ± 10.4 


1099.80 

0.448 

148.0 ± 13.7 

103.1 ± 5.4 

92.1 ± 6.8 

95.2 ± 3.9 


1112.59 

-0.024 

PF 

PF 

PF 

PF 


1115.91 

-0.081 

110.9 ± 9.2 

96.7 ± 6.3 

67.8 ± 3.7 

83.8 ± 4.0 

H 2 j - 4 . 

935.969 

1.264 


75.2 ± 6.4 




979.808 

1.095 


76.7 ± 7.1 

55.2 ± 3.8 



994.234 

1.134 


77.3 ± 8.5 

53.2 ± 6.2 



999.272 

1.217 


88.8 ± 8.1 

54.2 ± 3.3 

59.4 ± 4.9 


1017.39 

1.002 

PF 

PF 

PF 

PF 


1032.35 

1.247 

100.7 ± 6.4 

89.0 ± 6.1 

60.9 ± 3.1 

66.2 ± 3.0 


1044.55 

1.206 

91.9 ± 4.3 

86.7 ± 4.2 


72.9 ± 3.6 


1047.56 

1.062 

93.9 ± 6.5 

79.9 ± 4.4 

54.6 ± 2.6 

64.1 ± 3.6 


1057.39 

1.135 

97.0 ± 4.5 

72.8 ± 4.2 

59.9 ± 4.1 

65.8 ± 3.2 


1060.59 

1.019 

78.3 ± 5.3 

81.2 ± 3.9 

57.1 ± 2.8 

63.0 ± 3.4 


1074.32 

0.923 

91.1 ± 5.6 

71.2 ± 3.7 

51.5 ± 3.3 

57.7 ± 5.2 


1085.15 

0.817 


70.5 ± 8.8 

45.3 ± 4.0 



1088.80 

0.752 

78.8 ± 5.4 

71.6 ± 4.0 

46.7 ± 3.1 

58.4 ± 2.8 


1100.17 

0.498 

62.5 ± 4.2 

73.2 ± 5.5 

43.6 ± 3.7 



1104.09 

0.461 


64.5 ± 3.8 














Table 3 


Continued 


Species 

A (A) 

log (/A) 

HD 73882 WA(mA) 

HD 192639 WA(mA) 

HD 206267 WA(mA) 

HD 207538 WA(mA) 


1116.03 

-0.060 

34.1 ± 4.5 

44.7 ± 4.8 

17.6 ± 2.9 

28.3 ± 3.4 


1120.26 

-0.069 

30.5 ± 3.0 

45.4 ± 5.9 

20.2 ± 2.5 

33.3 ± 3.1 

H2 J - 5. 

942.691 

0.765 


37.5 ± 8.7 




974.889 

1.138 


68.4 ± 6.3 

38.8 ± 3.8 



996.129 

1.102 


PF 

45.4 ± 6.3 

40.2 ± 8.0 


997.644 

1.110 


62.3 ± 6.4 

44.5 ± 4.8 

33.3 ± 10.1 


1006.34 

0.940 

PF 

PF 

PF 

PF 


1017.01 

1.060 

88.3 ± 9.0 

67.6 ± 4.5 

PF 

PF 


1017.84 

1.384 

117.6 ± 7.5 

82.0 ± 3.8 

51.4 ± 3.1 

56.4 ± 5.0 


1040.06 

1.074 

91.2 ± 5.4 

64.9 ± 3.4 

37.1 ± 2.7 

33.7 ± 3.6 


1052.50 

1.068 

80.8 ± 6.2 

64.8 ± 3.2 

38.4 ± 3.4 

30.7 ± 3.0 


1061.70 

1.126 

87.8 ± 12.0 

76.8 ± 3.4 

39.1 ± 2.5 

45.4 ± 4.3 


1065.60 

1.026 

93.5 ± 9.2 

67.9 ± 3.6 

33.0 ± 3.2 

27.8 ± 3.2 


1075.25 

0.999 

94.9 ± 6.9 

70.6 ± 2.9 

36.6 ± 2.8 

31.6 ± 3.5 


1089.52 

0.796 

68.9 ± 4.3 

58.1 ± 4.6 

31.2 ± 3.4 



1104.55 

0.471 


30.9 ± 3.0 




1109.32 

0.467 

47.3 ± 4.1 

31.4 ± 3.1 

16.6 ± 2.0 

12.4 ± 2.2 


1120.41 

-0.095 

19.3 ± 3.3 

14.7 ± 2.7 

5.0 ± 0.9 

PF 

H2 J - 6. 

959.163 

1.426 


16.4 ± 4.2 




977.733 

1.538 


20.1 ± 4.2 




998.340 

1.564 


22.6 ± 4.2 




1019.02 

1.318 

31.4 ± 4.9 

14.0 ± 2.4 




1021.22 

1.388 

22.7 ± 6.9 

20.2 ± 3.6 




1041.74 

1.243 

26.1 ± 3.9 

11.8 ± 2.9 




1045.81 

1.076 

16.1 ± 5.1 

15.4 ± 3.5 




1058.32 

1.074 

18.0 ± 3.8 




H2 J - 7. 

1060.04 

1.193 

20.0 ± 3.9 





1073.00 

1.110 

14.2 ± 3.2 





Note. 


ines blended but used for profile fitting are noted “PF”. Errors are 1 a. 
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Table 4. COG determined b values for different combinations of EqW towards HD 192639 


Level 

All the EqW 

Half of the EqW 

Second half of the EqW 

H2 J = 2 

3-8 i 2!3 

4.3 ± 11 

o o _(_ 2.1 
^ 2.6 

H2 J = 3 

4.8 ± 0-.7 

5.0 ± li 

4-5 ± 1:1 

H2 J = 4 

5.6 ± 

5.4 ± O f 

6.0 ± O.T 

H2 J = 5 

8.0 ± 1:° 

7.8 ± 2-2 

7-9 ± 1:^ 


Note. 


b values are in kms Errors are 2 a. 







Table 5. 


Log column densities and b values obtained with PF and COG for the four lines 

of sight 


J level 


HD 73882 


HD 192639 


HD 206267 


HD 207538 


cm~^ / kms~ 

1 

COG 

PF 


COG 

PF 


COG 

PF 


COG 

PF 


V(H2) 

— J = 

0 

20.99 ± 

0.08 

0.08 


20.28 ± 

0.05 

0.05 


20.64 ± 

0.03 

0.03 


20.64 ± 

0.07 

0.07 


V(H2) 

— J = 

1 

20.50 ± 

0.07 

0.07 


20.48 ± 

0.05 

0.05 


20.45 ± 

0.05 

0.05 


20.58 ± 

0.05 

0.05 


V(H2) 

— J = 

2 

18.76 ± 

19.04 ± 

0.10 

0.09 

18.66 ± g:0i 

18.75 ± 

0.10 

0.05 

18.39 ± g:g| 

18.49 ± 

0.05 

0.04 

18.63 it g;Qy 

18.75 it 

0.05 

0.05 

V(H2) 

— J = 

3 

18.22 ± 0:10 

18.49 ± 

0.15 

0.04 

17.84 ± g;!! 

18.00 ± 

0.15 

0.10 

17.91 ± g:gi 

17.95 ± 

1.05 

0.05 

17.65 ± o-ll 

17.90 it 

0.10 

0.15 

V(H2) 

— J = 

4 

15.66 ± 0.10 

15.80 ± 

0.10 

0.16 

16.01 ± o-.u 

15.90 ± 

0.10 

0.10 

15.46 ± Q:j^2 

15.28 ± 

0.17 

0.10 

15.77 ± g:i| 

15.70 it 

0.10 

0.06 

V(H2) 

— J = 

5 

15.34 ± 0.11 

15.34 ± 

0.10 

0.11 

15.18 ± 0-^7 

15.15 ± 

0.08 

0.07 

14.83 ± g:ig 

14.75 ± 

0.15 

0.10 

14.59 it g;Qg 

14.64 it 

0.16 

0.15 

V(H2) 

— J = 

6 

14.28 ± 0.11 

14.08 ± 

0.26 

0.08 

13.98 ± g;!! 

14.00 ± 

0.15 

0.10 


13.84 ± 

0.15 

7 


13.04 it 

0.60 

7 

V(H2) 

— J = 

7 

14.16 ± 0.17 

14.08 ± 

0.37 

0.23 


13.49 ± 

0.35 

7 


13.59 ± 

6.31 

7 


13.28 it 

0.56 

7 


5(H2) 

— J = 

2 


2.47 it f g 

3-8 ± 2.3 

4 KQ _|_ 0.40 

4.oy m Q 7 Q 

3.8 ± 

9 40 _L 0.20 

ziz Q 4 Q 

3.3 ± i:i 

3.97 it 

0.20 

0.20 

dH 2 ) 

— J = 

3 

5-7 ± fi 

4.39 it 

4.8 ± g:g 

5.14 ± g:ig 

3-5 ± g:i 

3.10 ± g:fg 

4.5 ± g:i 

4.41 it 

0.20 

0.20 

5(H2) 

— J = 

4 

7-7 ± g:^ 

7.19 ± g;| 

5-6 ± g:g 

5.89 ± g:ig 

4.7 ± g:t 

5.02 ± g:li 

4.8 ± g:« 

4.78 it 

0.20 

0.30 

b(H 2 ) 

— J = 

5 

-H 

o 

11.09 it fg 

8.0 it i:g 

8.34 ± g:tg 

5.3 ± I;! 

6.50 ± 0:tg 

10.8 it I;? 

7.22 it 

2.40 

1.60 


Note. 


J = 0 Sz 1 column densities from Rachford et al. (2002). Errors are 2 a. 
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Table 6. Gas Densities and H 2 Formation rates — b due to H 2 -formation pumping 



n 


nR 


R 

Lines of Sight 

[cm“^] 


[s-i] 

[cm^ s 1] 

HD 73882 

5± f 

9± 

8 X 10-15 

1.8 ± 

1.2 X 10-15 

HD 192639 

HD 206267 

13± i 
25 ± 

4± 

6± 

5 X 10-15 

1 X 10-15 

3.1 ± 

2.4 ± 

X 10-15 

X 10-15 

HD 207538 

13 ± 

3± 

2 X 10-15 

2.3 ± 

1.9 X 10-15 


Note. — Densities and formation rates come from our mea¬ 
surements reported in Table 5 and model results plotted in 
Figure 7. 







Table 7. 


Temperatures and Densities — b due to thermal scattering 




n~'ex 

-^J=b-7 

n 

logm{P/k) 

Lines of Sight 

(K) 

(K) 

(cm“6) 

(log(K cm“^)) 

HD 73882 

9500 ± 600 

670 ± Z 

31 ± 

4 3+ 0'^ 

0.3 

HD 192639 

8200 ± 

500 ± po 

13 ± f 

3.8+ 

HD 206267 

7000 ± 

680 ± fo 

49 ± f 

4.5+ V 

HD 207538 

8500 ± ?»»„ 

610 ± fo 

27 ± f 

4.2+ r 


Note. — is the temperature obtained from the formula b = 
\/2kT/m. is the excitation temperature from the rotational 

levels 5 and 7. 







